Selection-rule blockade and rectification in quantum heat transport 
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We introduce a new thermal transport phenomenon, a unidirectional selection-rule blockade, 
and show how it produces unprecedented rectification of bosonic heat flow through molecular or 
mesoscopic quantum systems. Rectification arises from the quantization of energy levels of the con- 
duction element and selection rules of reservoir coupling operators. The simplest system exhibiting 
the selection-rule blockade is an appropriately coupled three-level system, providing a candidate for 
a high-performance heat diode. We present an analytical treatment of the transport problem and 
discuss how the phenomenon generalizes to multilevel systems. 
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Heat conduction in nanoscale structures has become an 
active field of research enjoying constantly increasing at- 
tention. Electronic properties have been studied in great 
detail in the last three decades and are better under- 
stood than thermal properties. The main reason for this 
has been experimental challenges to measure and con- 
trol thermal properties accurately. Recently the field has 
seen breakthroughs such as measurements of quantized 
heat transport [TJ [2] and realization of hybrid structures 
probing atomic-level heat transport properties [5] . 

Experimental developments have increased theoretical 
interest to explore fundamental limits of thermal phe- 
nomena and devices based on their applications. Ther- 
mal rectification, i.e. asymmetry of heat current when 
the temperature bias is inverted, has been actively stud- 
ied in this context. Rectification has been observed so far 
in two experiments |U |S] and has a significant applica- 
tion potential. There exists various theoretical proposals 
of how to realize rectification in phononic [5J [7], pho- 
tonic [3] and other hybrid structures !! . However, in 
most proposals rectification is modest and limits of per- 
formance are unknown. The purpose of this paper is to 
introduce a novel heat transport phenomenon, an asym- 
metric selection-rule blockade, enabling an unparalleled 
rectification performance. The phenomenon provides an 
example of an intricate interplay of quantum-mechanical 
and thermal properties in nanoscale structures. 

The simplest system exhibiting the selection-rule 
blockade is a three-level system with a strongly uneven 
energy-level separations. A crucial ingredient is that the 
two baths are coupled to the system so that one can only 
induce transitions between the close-lying states and the 
other can effectively induce only the two larger transi- 
tions. The baths can exchange energy only when the 
bath coupling far-apart states has sufficiently high tem- 
perature to create excitations. In suitable temperature 
regime the forward biasing leads to sequential heat flow 
while the reverse bias current, due to higher-order pro- 
cesses, is very strongly suppressed. This one- directional 
suppression of heat flow is the defining property of the 
selection- rule blockade. 



The paper is organized as follows. First we introduce 
the three-level model and show analytically how rectifi- 
cation arises from the selection rules. Then we estimate 
the magnitude of the leakage current that sets the limits 
of performance of rectification and discuss how the se- 
lection rules can be realized without special symmetries. 
We conclude by outlining how the selection-rule block- 
ade generalizes to multilevel systems and summarize our 
results. 




FIG. 1: Studied model consists of two heat baths at different 
temperatures coupled through a three-level system. The blue 
arrow represents transitions induced by the left reservoir and 
the red arrows correspond to Golden-Rule transitions induced 
by the right reservoir. The dashed arrow represents higher- 
order transitions via the virtual state by the right reservoir. 

In the following we consider thermal conduction in the 
setup depicted in Fig. [l] The system consists of three 
parts, two reservoirs and a three-level system that me- 
diates heat. We call a configuration forward biased if 
Tr > XL and reverse biased in the opposite case. Intu- 
itively one expects that if temperatures Tl/r are much 
smaller than the allowed transition energies, heat flow is 
effectively blocked. This is confirmed below as we show 
how the intrinsic asymmetry in the strength of the block- 
ade gives rise to strong rectification. The Hamiltonian of 
the system is 

H = H L + H R + H C + H LC + H RC (1) 

where -Hl/r characterize the reservoirs and the central 
Hamiltonian is 

ff C = eo|0)<0| + e 1 |l)(l| + e 2 |2)<2|. (2) 
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The level separation of the central part E\ = t\ — eo- 
E2 = £2 — ei is assumed to be unequal E1/E2 <C 1. 
The reservoirs couple to the central system in a special 
manner 



ffLC=*L(|0)(l| + |l)<0|) 

Hrc = ^r(|1>(2| + |2>(1|- 



|0)(2| + |2)(0|), (3) 



where ^l/r ar e Hcrmitian operators in the reservoir 
Hilbert spaces. The specific form of the couplings im- 
ply that the left reservoir couples only the groundstate 
and the first excited state in contrast to the right reser- 
voir which couple all the energy levels. Applying the 
Fermi Golden Rule one can calculate the transition rates 
between the energy levels 

Ti^o = Sx^ij/h 2 , r 2 ^i = Sx R (w 2 )/ft 2 , 



(4) 



where Sx L/t » = dte*-*(X L/R (t)X L/R (0)) and 

u)\/2 = E1/2/K. The correlation functions are evaluated 
in the absence of the couplings at the respective reser- 
voir temperatures. Information about Hi, and i? R as 
well as information about the structure of the reservoirs 
is encoded in the noise Sx L / R (w) so our discussion is in- 
dependent of the precise nature of the reservoirs so far. It 
should be noted that even though the right reservoir can 
induce direct transitions between the ground state and 
the first excited state it happens only in higher-orders. 
In the lowest order the allowed transitions are separate in 
the different reservoirs. The corresponding inverse tran- 
sitions are obtained by inverting the sign of the frequency 
argument, for example r ^i = Sx L ( — ^l)/^ 2 - Noise is 
related to the retarded X-correlation function and the 
Bose-Einstein distribution n(u>) through the Fluctuation- 
Dissipation Theorem 



where Al/r(w) = — 2 Im(X L / R X L / R ) r (tL>) are spectral 
densities of the reservoirs defined through the Fourier 
transform of the retarded function (X L / R (i)X L / R (0)) r = 
-iO{t)([X L/R (t),X L/R (0)]} [ID]. Once the transition 
rates are known, the steady-state occupation probabil- 
ities follow from the detailed-balance equations 



0. 



(6) 



Employing the rate-equation formulation |10j . heat cur- 
rent can be calculated, say, between the left reservoir and 
the central system yielding 



J = Ei [P(l)ri^ - P(0)IW] 



(7) 



The probabilities required to evaluate Q can be solved 
in terms of the transition rates as 

p(o) = (r 1 ^ r 2 ^ + r 1 ^ 2 r 2 ^ + r^cTW) /c 
p(i) = (iwr 2 ^o + iWiw + r ^ 2 rw) /c, (8) 



with 

c = Tioi^o + r 12 r 20 + r 10 r 21 + r i (r 20 + r 12 + r 21 ) 
+ r 02 (r 21 + r 10 + r 12 ), (9) 

where we have introduced notation Ti—>j = Tij . Inserting 
the probabilities and the transition rates in Eq. ^ we 
obtain an explicit expression for the current 



J = E t 



roiFo2r 2 i 



c 



(10) 



Result (10) indicate that the thermal window [eP hEl 



e^ nEl ) determined by the reservoir temperatures and the 
smaller energy separation E\ plays an important role in 
the transport process. In the limit /cb?r ^ E2 current 




FIG. 2: Rectification efficiency as a function of the higher 
operation temperature. The different curves correspond to 
values T = Si, and E2/E1 = 5 (red), E 2 /E 1 = 7 (blue), 
E2/E1 = 9, (green) and £ 2 /-Ei = 11 (black). Inset illus- 
trates the corresponding forward bias (solid lines) and re- 
verse bias currents (dashed lines, almost zero) in the units 
of Jo = EiA /h 2 . 



(5) ( 10 1 takes a simple form 



1 



A r (uj 2 ) 



(l + e&Ei) 



(11) 



Result (111 shows that heat current is exponentially sup- 
pressed in the considered limit. This is an indication that 
the system acts as a high-performance thermal rectifier 
when the lower operation temperature is much smaller 
than E2 while the higher temperature is comparable to it. 
To make rectification more explicit, consider a situation 
where the hot reservoir is at temperature /.'bHi ~ Eo and 
the cold one is at k^T c ~ E\. Furthermore, to simplify 
expressions, let us assume that the effective couplings are 
of the same order A\,{u}\) w j4r,(w 2 ) w Ah(wi+w 2 ) = Aq. 



Starting from the expression ( 10 1 one can estimate that 



the fraction of the bias-inverted heat currents is 



J T R =T h 



-E 2 /Ei 



-T h /T c 



(12) 
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where c is a numerical factor of the order of unity. Thus tribution becomes 



Eq. ( 12 1 shows that in the studied temperature regime 



the fraction of forward and reverse bias currents is expo- 



nentially small in T^/T c . Results ( 10 >-( 12 ) illustrate the 
essence of the selection rule blockade and confirm quanti- 
tatively the intuitive physical picture. More comprehen- 
sive picture of transport properties can be obtained by 



plotting current (10 1 and the rectification ratio at differ- 



ent forward and reverse bias values, see Fig. [2] Reverse 
bias current is rapidly suppressed as the ratio E 2 /Ei 
is increased and effectively vanishes in the temperature 
window fceT c < E\, fceTh > 3£a for E 2 jE\ > 3. Impor- 
tantly forward bias current is finite and increasing in the 
region of optimal rectification. 



In the derivation of Eq. (10 1 we only took into ac 



count the Golden-Rule transitions. Deep in the block- 
aded regime fee^R <SC E 2 the lowest-order current van- 



ishes exponentially as shown in Eq. (Ill and one is nat 



urally lead to consider higher-order transitions. The sit- 
uation is similar to electronic transport in the Coulomb 
blockade model where charging effects suppress conduc- 
tance (with the difference that the selection-rule blockade 
is one directional). At low temperatures and biases elec- 
tric current through a small island vanishes exponentially 
in the lowest order in the perturbation theory. However, 
in the next order one obtains cotunncling processes where 
electrons are transferred via virtual intermediate states 
partly lifting the blockade [TTj . In our model the leading 
higher-order processes contributing to the reverse leak- 
age current in the blockaded regime correspond to the 
right reservoir-induced transitions via the highest-energy 
state as depicted in Fig.[T] These rates can be calculated, 
for example, by the standard T-matrix expansion to the 
second order [10]. Let us assume that the right reservoir 
consists of noninteracting bosonic modes and that the 
coupling operator is of the form Ar = Yljen c j{bj + bj) 
where bj , bj are canonical bosonic operators of the mode 
ujj. In the low-temperature limit Tr — > the relevant 
second-order process contributing to reverse bias current 
is the decay of the first excited state to the ground state 
by emitting two bosons to the right reservoir. This pro- 
cess gives rise to an additional rate 



r (2) 

1 1->0 



2tt 



du>An(u>)Aji(LJi — lu)x 



1 



1 



Ei — huj E 2 — Ei + huo 
8tt r ut 



KE\ jo 



a\jjAn(uj)An(ui —u>). 



(13) 



The rate is proportional to a typical second-order energy 
denominator ~ 1 / E 2 suppressing the transition. As- 
suming that the spectral density has a power-law form 
Aa(w) cx uj d (d — 1 for an ohmic bath), the leading con- 



,(2) 
l->0 



in (hT 2 ^Q\T B =oY 



1 + d 



E 2 



2d 



(14) 



The second and the third factor are much smaller than 



unity so ( 14 1 is small compared to the level separation 
u>i. However, the rate is only algebraically suppressed by 
(Ei/E 2 ) 2d . At finite temperature Tr ~ E\ the relevant 
power law exhibits a crossover to (fcsTfl/^) . Now we 
can estimate the leakage current under the blockade by 



J r = -£^(1)1^0 = -Ei 



r (2) 

1 I->0 

r (2) 

1 + e A^i + 

1 l->0 



(15) 



Result ( 15 1 shows that the true rectification efficiency is 
algebraically, not exponentially, suppressed. It is clear 
that the leakage current remain finite as long as there is 
thermal coupling between the two reservoirs and that a 
good rectifier is characterized by its ability to suppress 
the parasitic reverse bias processes. Nevertheless, cur- 



rent (15 1 can be made arbitrarily small by decreasing 



max{E 1 /E 2 ,k B T K /E 2 }. 

Usually selection rules are associated with symmetries 
of the Hamiltonian leading to vanishing matrix elements 
of the perturbation operator between unperturbed states. 
In applications finding a candidate system with suitable 
symmetries and couplings is a nontrivial task. However, 
since the transitions of the different reservoirs are well 
separated in energy, there exists an alternative route to 
obtain the selection rules without invoking symmetries of 
the Hamiltonian. Golden-Rule rates consists of matrix el- 
ements of perturbation operators and a summation over 
a relevant phase space. Both factors can effectively im- 
pose selection rules in the system. Typically phase-space 
selection rules arise from the fact that the reservoir spec- 
tral density is nonvanishing only in a restricted frequency 
window outside which the reservoir cannot induce tran- 
sitions. In solid-state applications suitable reservoirs can 
be realized by dissipative vibrational modes of phononic 
or photonic nature which filter out frequencies far from 
the resonance. Assuming that the center element cou- 
ples linearly to the modes, the reservoir parts of the cou- 
plings |i| take the form X L / R = c L / R (6^ r +& l /r), where 

^l/R' bh/R ar e the creation and annihilation operators of 
the reservoir modes. The retarded Green's function of a 
vibrational mode is G r {uo) = ((fift) -1 + ^ r ( a; )) _1 j where 
(ffo) -1 = ( w2 — t^oV^o is the inverse free Green's func- 
tion, u>q is the frequency of the mode and S r (w) is the re- 
tarded self-energy due to dissipation. Dissipation can be 
realized by coupling the mode to a bosonic bath, in which 
case the spectral densities Al/r(w) = — 2gf j ^ R lmG r (uj) 
of the reservoirs can be solved exactly [12] yielding 



c l /r( 2 ^l/r)X/r( 



^ (w 2_ w 2 /r)2 + (A B /rMwl/r)2 ' 



(16) 
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where wl/r are the resonant frequencies of the reservoirs 
and Aj^y R (a;) are the spectral densities responsible for dis- 
sipation of the reservoir modes. As illustrated in Fig. [3] 
(left), if the left reservoir has a resonance frequency at 
ujt, = u>i and the right reservoir at ojr = ^2 + wi/2 with 
respective widths |15j that are much smaller than u>2, 
the spectral densities have no overlap. Thus the reser- 
voirs can effectively induce transitions only as indicated 
in Fig. [T] If the reservoirs are electromagnetic in nature 
the desired spectral densities (16 1 are achieved by cou- 



pling the center clement to dissipative LC-circuits whose 
electromagnetic fluctuations are restricted to a narrow 
band around the resonance [13 . In phononic systems 
the same effect could be realized by coupling large baths 
to the central clement through small vibrating bridges. 




FIG. 3: Left: Spectral density of the left (right) reservoir 
is completely suppressed at frequencies u»wi (uj <§C lo 2 ), 
effectively imposing the desired selections rules. Right: Low- 
lying spectrum of the Jaynes-Cummings model and the al- 
lowed transitions induced by the left (blue arrows) and the 
right (red arrows) reservoir. 



The selection-rule blockade is a not restricted to the 
studied three-level system. Basic requirements are that 
level separations of a system can be divided to small 
and large intervals coupled to separate baths, and that 
combined effect of two large transitions corresponds to 
a small one as in the tree-level model. A prominent ex- 
ample fulfilling the requirements is a two-level system 
coupled to a harmonic oscillator described by a resonant 
Jaynes-Cummings- type Hamiltonian |14j 



H 



jc 



Ml [ a ' a - 



Ml Kg + 

— cr z + —{a<+a)a x , (17) 



where the interaction is assumed to be small g <gc fi. Ex- 



cited states of the Hamiltonian (17) consists of doublets 
\E n ) — \n, I) ±\n— 1, t) where the energy separation of 
centers of adjacent doublets is approximately Eq = Ml 
and the separation within a doublet is E n = ^nhg 7 see 
Fig. [3] (right). Since for the low-lying doublets energy 
splittings are small E n <C E , the selection-rule blockade 
can be established if the left reservoir couples levels only 
within doublets and the right reservoir couples states in 
different doublet. This requirement is satisfied in the 
lowest order by H RC — X R (a^ + a) and #lc = X L a z , 
yielding immediately the desired selection rules. Anal- 
ogous to the three-level model, if the temperatures of 



the reservoirs are A;bTl ~ %j &bIr ~ Ml heat flow is 
efficient, while in the reverse biased case /cbTl ~ f&l, 
fcB^R ~ % excitations induced by the right reservoir are 
forbidden and there exists only a weak leakage current 
due to higher-order processes. 

In conclusion, we introduced a fundamental heat trans- 
port phenomenon, the selection-rule blockade, based on 
the quantum nature of the central system and selection 
rules of the bath coupling operators. The phenomenon 
enables a high rectification efficiency in molecular and 
mesoscopic structures making it promising for future de- 
vice applications. The simplest system exhibiting an 
asymmetric blockade is a three-level model which we ana- 
lyzed in detail. We proposed a scheme to impose required 
selection rules without special symmetry properties of the 
Hamiltonian and discussed how the selection-rule block- 
ade can be generalized to multilevel systems. 
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